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Abstract. 



We use the cosmic microwave background temperature anisotropy to place limits on large-scale 
magnetic fields in an inhomogeneous (perturbed Friedmann) universe. If no assumptions are made 
about the spacetime geometry, only a weak limit can be deduced directly from the CMB. In the special 
case where spatial inhomogeneity is neglected to first order, the upper limit is much stronger, i.e. a 
few xlO" 9 G. 



1. Introduction 

Magnetic fields have been observed in the universe on a wide range of scales. Fields with strengths of a 
few /iG are prolific in galaxies and galaxy clusters, extending well beyond the core regions of the latter, 
and have also been detected in high redshift Lyman-a objects. Magnetic fields in extragalactic structures 
are detected mainly via radio polarisation studies, X-ray emission and Faraday rotation measurements 
(see [1] for a comprehensive review). Magnetic fields in galaxies and galaxy clusters appear to be the 
result of the nonlinear amplification of weak seed fields, mainly via the galactic dynamo. However, the 
detection of ordered magnetic fields in high redshift objects (with z > 2) poses a stiff challenge to the 
dynamo mechanism. As yet, there is no direct evidence of magnetic field presence on cosmological scales, 
corresponding to a significant fraction of the Hubble length. Clearly, any such field could not arise through 
structure formation physics, but it would have to be the remnant of a primordial field, redshifting with 
expansion: 



where Bq is the current field strength. 

The strength of primordial, cosmological magnetic fields is limited by observed helium abundances 
and by the near-isotropy of the cosmic microwave background (CMB) (see [2] for recent reviews). Any 
magnetic field present at the time of cosmological nucleosynthesis inevitably affects the abundance of 
primordial Helium, since it provides an additional form of relativistic energy density. This, in turn, 
increases the expansion rate of the universe with the effect that the neutron-proton freeze out of weak 
interactions occurs at a higher temperature. The result is an increase in the synthesised abundance of 
primordial Helium. Hence, Helium-4 observations (extrapolated to zero metalicity) provide an upper 
limit of ~ 10~ 7 G, in today's values, on any primordial magnetic field present at nucleosynthesis [3]. 
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Stronger limits are imposed from the observed high isotropy of the CMB photons. The COBE 
data place an upper bound on a homogeneous magnetic field present at the time of last scattering. In 
a recent analysis of a particular class of spatially homogeneous Bianchi universes, an upper bound of 
B Z 10~ 9 G was obtained [4]. Here we generalize previous work to the case of inhomogeneous fields in 
an inhomogeneous almost-Friedmann universe. We also generalize the limits found in [4] by weakening 
some of their assumptions. It turns out that, if we do not assume a spatially homogeneous geometry, the 
limits imposed directly by CMB data on super-Hubble magnetic fields in an inhomogeneous universe are 
much weaker, Bo S 10~ 6 G. A similar situation arises when considering the limits placed on the shear by 
CMB anisotropies, as pointed out in [5]. 

Recently it was proven [6] under quite general circumstances that a magnetic field is prohibited in 
spacetimes where exactly isotropic radiation is also present. Taking this as our starting point, small 
anisotropy allows for a weak magnetic field. We use the 1+3-covariant analysis of CMB temperature 
anisotropies [7, 8, 9] and of magnetic fields [10], in order to derive limits on large-scale fields as a function 
of coherence scale. Following the approach of [8, 11], we use the radiation multipoles to derive limits 
which are model-independent, in the sense that they do not rely on assumptions about the (perturbed 
Friedmann) spacetime geometry. 

In Sec. 2 wc outline the general formalism for imposing limits on large-scale magnetic fields from 
observed CMB temperature anisotropics. In Sec. 3, we give our main results, which follow from a 
refinement of the method in [11]. For convenience, we omit most of the calculational details, and give 
the key equations in Appendices A and B. We use units such that 8ttG = 1 = c. Our notation follows 
that of [12]. In particular, X a ... b = u c V c X a ... b and D a is the covariant derivative in the rest space, i.e., 
D c X a ... b = h c d ha ■ ■ -h/V 'dX e ...f, where h ab = g a b + u a u b is the projection tensor. Angled brackets on 
indices denote the projected, symmetric and trace free (PSTF) part. The 3-divergence and 3-curl of 
PSTF tensors are divX a = B b X ab , ■ ■ •, and cur\X a = e abc D b X c , curl X ab = e cd(a D c X b) d , 

2. CMB anisotropy induced by large-scale magnetic fields 

In the 1+3-covariant analysis of CMB anisotropy [8, 9], a physical choice of 4- velocity u a is made, usually 
the 4- velocity of cold dark matter (CDM), and all perturbative quantities are then covariant vectors or 
tensors in the rest-space of u a , with direct geometrical or physical meaning. The fractional temperature 
fluctuation is expanded in covariant multipoles TA e (Ag = a\ ■ ■ ■ a?), which are PSTF tensors. These are 
limited directly by observations: 



which we use here. The first moment r a is the dipole, which is usually attributed to our peculiar motion 
relative to the CMB frame. We assume this motion is corrected for by setting r a = = ei for the bulk 
of the paper. However, it is possible that a residual dipole of cosmological origin exists (it would be 
frequency dependent, and thus could not be set to zero by a Lorentz boost), and we include it in our 
calculations for generality. 

In addition to the observed bounds on the t A( , we need bounds on their temporal and spatial 
gradients, in order to find limits on geometrical and physical quantities that characterize the spacetime. 
Wc define the expansion-normalized, dimcnsionlcss e-quantitics [8, 11] 
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COBE data leads to the values [13] 

e 2 w 1.1(±0.8) x 10 
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where H is the background Hubble rate. 
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Following [11], we can find upper bounds on all perturbative quantities in terms of the e e , e£, e\, 
etc. However, the derivative bounds e£, e\ ■ ■ ■ are not directly measurable, as we are unable to move 
a cosmological time- or space-separation from our current spacetime event. Here we make the simple 
assumption that the time- and space- variations of the multipoles are governed respectively by the Hubble 
rate and the physical scale A of the perturbation, i.e., 

\t Ae \ ~ H\t A( \, \D a TA t \ ~ j\ta ( \- (6) 
This assumption implies 

4 ~ g 6 ^ e T ~ g e <?, ' ' ' ! ( 7 ) 

cj ~ ^Q, 4 f ~ g^e/. • • • where P = jpr- ( 8 ) 

The dimensionless parameter gives the coherence scale as a fraction of the Hubble length, with £ 0(1) 
since we are considering large scales. 

The magnetic field is 'frozen' into the baryonic fluid, which may be treated as an infinitely conducting 
medium. Here, we neglect the peculiar velocity of CDM relative to the baryons. The physical justification 
for doing so comes from the fact that we address the linear regime and consider large scales, where the 
velocity difference between the two components is expected to be minimal. On these grounds, we choose 
u a to be the 4- velocity of CDM-baryon fluid with total density pm- The kinematics of the pressure-free 
matter are characterized by the volume expansion 6, rotation co a , acceleration A a and shear distortion 
o- a b of u a . Note that, even in the absence of pressure gradients, the flow lines are generally non-geodesic 
(i.e. A a 0) due to the magnetic field presence. Here, however, we will assume an effectively force-free 
field (i.e. e a bcB b cvs\B c = 0), and ignore the acceleration to first order. This is a reasonable approximation, 
given that the field is too weak to affect the motion of the baryonic matter. It will also allow us to focus 
upon the purely anisotropic magnetic effects. 

The energy-momentum tensor of the magnetized dust is 

Tab = (pM + PB)u a U b + p B h a b + ^ab , (9) 

where p& = B 2 /8ir and pb — Pb/3 are the magnetic energy density and isotropic pressure respectively. 
Also, H a b = —B( a Bty/&-K is the symmetric and trace free tensor that conveys the anisotropic magnetic 
effects. The radiation energy-momentum tensor is 

Tab = pU a U b + ^fih ab + 2lt( a <7 b ) + 7T a6 , (10) 

where p, q a and TT a b are the photon energy density, momentum density and anisotropic stress. These are 
directly related to the temperature anisotropy multipoles by [8] 

= »PT a , TTab = T7V T ab- (H) 

o 15 

The photon energy momentum tensor involves only the first two multipoles, but we will require also the 
octupolc 

g 

iabc = T^P-Tabc, (12) 

which appears in the evolution equation for Tt a b, Eq. (A. 7). 

The field equations G a b = T ab + T ab , the Ricci identities and the Bianchi identities may be split into 
a set of evolution (along u a ) and constraint equations. The evolution of the magnetic field is determined 
by Maxwell's equations. The reader is referred to Appendix A for the necessary equations. 



3. The limits 



We first present the CMB limits on inhomogeneous magnetic fields as a function of coherence scale, A; 
then we discuss the homogeneous case. 
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3.1. Inhomogeneous universe 

Our procedure to find constraints on the magnetic field strength p Bl or equivalently |II afc |, is a 
generalization of the non-magnetized analysis in [11]. Briefly, we manipulate the field equations to 
express Tl a b in terms only of the radiation quantities fx, q a , n a b, £ a bc- This is facilitated by the appearance 
of the shear in Eq. (A. 7), which, in the absence of acceleration, is the only coupling of the radiation to 
the first-order kincmatical quantities. The main aspects of this calculation are in Appendix B, and the 
key result is Eq. (B.8). 

Neglecting the dipole moment and the energy density of the radiation Or, and restoring units, 
Eq. (B.8) gives 

3\ 3/4 ctf f (*_<, , 45_ t \ , n ( , 9_ t 



Bo < max(S) =^-j ~y= [5e* + — e\ j + A l^e* + -e 3 

18 t 135 t 81 t * 81 m l 1/2 

7 3 14 3 14 3 14 3 J ' [ ' 

where the e's are evaluated at the current time. The function max(B) gives upper limits on large-scale 
magnetic fields, coherent on a given scale A, imposed by CMB temperature anisotropics. This upper 
limit is given directly in terms of CMB multipoles and their derivatives, and is model-independent, i.e., 
no assumptions have been made about the spacetime geometry. 

For a numerical estimate, we need to use the simple assumptions of Eqs. (7) and (8), in order to 



nmx(^,= |'|) -9H<2 



+£3 



we 


find, 


3/4 


cH 




Vg 


"1 








7 





5 15 1 1 

12 3 3 2 A 1(T 



M 
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> . (14) 



This is our main result. Using the limits in Eqs. (3) and (4), we can evaluate max(i3), which is plotted 
in Fig. 1. One of the main features of the plot is that uncertainty in £2 and 63 from COBE data produces 
a far greater uncertainty than the uncertainty in the cosmological parameters. 
On the largest scales, (3 — > 00, Eq. (13) simplifies to give 

3 / 4 „tt. ( k 1 =1 V2 



—<*)+. - (|) 9{|n„ + In A + !} e2 ~io-G. 



(15) 



3.2. Spatially homogeneous universe 

The upper limit on Bq on the largest scales in the general case of an inhomogeneous universe, as given by 
Eq. (15), is much weaker than the limit that can be imposed if we assume that the universe is spatially 
homogeneous to first order. Homogeneity implies that we can set to zero the e's that involve gradients, 
as we did in deriving Eq. (15). + However, it is not only the radiation multipoles that are homogeneous 
to first order, but the whole spacetime, leading to a Bianchi model. The special dynamics of Bianchi 
models then leads to a tighter constraint on B . A similar situation arises when deriving limits on the 
shear a a b [5]. 

It follows from [11] that the spatial 3-curvature vanishes to first order, 1Z a b = Tl = 0. In addition 
the shear becomes, from Eq. (A. 7), 

<Jab = -Tab- (16) 
+ Note that homogeneous radiation multipoles, i.e. D a TA e = 0, imply a Bianchi spacetime to first order [8]. 
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Figure 1. Maximum magnetic field strength in fiG on large scales as a function of coherence 
scale A (> 1000 Mpc) for various cosmological parameters {S2m^a} shown in the key, with h = 
-ffo/(100 km/s/Mpc) = 0.7. The concordance model {S!m^a} = {0.3,0.7} (central thick solid curve) 
is straddled on either side by the collective uncertainty in £2 and 63 (outer thin solid curves). An 
open and a closed model with extreme parameters are shown by the dashed curves, indicating how 
little the cosmological parameters affect the upper limit, and certainly much less than the observational 
uncertainty. The vertical line represents the Hubble distance, A = 3000/i — 1 Mpc. Because our model is 
only valid down to scales of the order of the Hubble distance, scales significantly below this should be 
interpreted with care. 



Thus the shear evolution equation becomes, using Eq. (A. 18), 

g 

Ilafc = Tab + Bf afc + — fXTab- (17) 

15 

The magnetic field acts as a forcing term for the quadrupole. The particular solution associated with 
this forcing term is 

g 

Ilab = TT^ab, (18) 

15 

with f a b = for the particular solution. The solution to the (non-magnetic) homogeneous part oscillates 
(at frequency ~ y/8fj,/15), while being suppressed with a damping scale of the Hubble time. Thus a 
conservative upper limit is given by Eq. (18); using Eq. (4) and Or ~ 2.5h~ 2 x 10~ 5 fi M) we find that 

B Q < 6.2±H x lO- 9 y/U^ G. (19) 

With fl M = 0.3, Eq. (19) gives 

B Q < 3A±H x 1(T 9 G. (20) 

This confirms the value found in [4], and is derived under slightly weaker assumptions; the spacetime 
is not chosen as a specific exact Bianchi model, but is homogeneous to first order, and turns out to be 
Bianchi I if we start with a flat Friedmann background. Furthermore, we include both the radiation 
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energy density and the cosmological constant. In the set of Bianchi models which admit a pure magnetic 
field (types I, II, III, IV Q , VII ), they are all of the same genericity; therefore we may consider this more 
general than [4] , where the geometry is assumed to be type VII/j . 

4. Conclusions 

For large-scale magnetic fields in an inhomogeneous almost-Friedmann universe, we have found upper 
limits on the field strength directly in terms of the CMB temperature multipoles and their derivatives, 
as given by Eqs. (13) and (14). On super-Hubble scales, this upper limit is very weak: 

B <2x 10~ 6 G 

for the concordance model, as shown by Eq. (15) and Fig. 1. 

When the almost-Friedmann universe is assumed to be homogeneous to first order, i.e. a Bianchi 
spacetime, the upper limit is much stricter, as given by Eq. (19). This generalizes the result of [4], by 
including a cosmological constant and removing initial assumptions of a choice of model. 

These limits have been derived in a covariant and gauge invariant way using the 1+3 formalism. 
A major feature of our approach is that our limits are largely model-independent, being derived from 
properties of the Einstcin-Boltzmann equations. Our main assumptions are imposed on the e-quantitics 
that bound the derivatives of radiation multipoles, which are in principle observable but in practice are 
not measurable. 

It is also possible to find limits on the inhomogeneity of the magnetic field, as given by the gradient 
of the magnetic energy density D a( o B - The result is given in Appendix B by Eq. (B.12). 
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Appendix A. The linearized equations 

We assume that CDM and baryons share the same 4- velocity, which coincides with that of the fundamental 
observers. Also, confining ourselves to times after last scattering, we may treat the magnetized dust and 
the radiation fluid as independently conserved entities (i.e., V b T a h = = V b T a i,). Then we arrive at the 
following linearized evolution equations for the magnetised dust (see [10]) 

PM + Qpu =0, (A.l) 

PM^a + ^DaPB + divn^O, (A.2) 

Pb + ^Opb =0, (A.3) 

iu + ^en ab =0, (A.4) 

4 

/i+ -O^i + divg =0, (A. 5) 

4 4 1 

4a + 3 e 9a+ jM« + -D a ^i + dlV7T a =0, (A.6) 

4 8 

TTab + gOTTab + J^^ab + 2D (a q h) + div£ afc = , (A.7) 



and 
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for the photons (see [12]). Note that in deriving Eq. (A. 2) we have used the linear relation e abc B b cur\ B c = 
^D a B 2 + D b II ab (recall that D a _B a = 0). To first order, the kinematic evolution is given by 

e + ie 2 + i(p M + 2M + 2p B )-div^-A =0, (A.8) 

2 11 
Oab + g©er ab + E ab - -ir ab - -U ab - B {a A b) = , (A. 9) 

LJ a + ^ecj a + ^cur\A a = 0, (A.10) 

2 

(?a - gD a 9 + div cr a - curl w a =0, (A. 11) 

divw =0, (A.12) 

where A a is given by Eq. (A. 2). Finally, the spacetime geometry is determined by [12] 

1 4 

E ab + QE ab - curl H ab + -(p M + ^)o- ab + 

+ ^9n ab - ien ah + ^B {a q b) = , (A.13) 

H ab + QH ab + curl E ab - ^curl ir ab - ^curl II ab = , (A. 14) 

H ab - curl cr afc - D (a w fc) = , (A. 15) 

div£ a + idiv7r a + idivn a - ^D (p M +M + Pb) + ^6<?a = 0, (A.16) 
1 4 

div H a + -curl q a - (p M + -p)uj a = , (A.17) 

T^-(ab) + gB(T a b - -7T ab - -j^-ab - E ab = , (A. 18) 

TZ-2(p M + fj, + p B ) + ^Q 2 -2A = 0, (A.19) 

where lZ ab and TZ are respectively the projected Ricci tensor and Ricci scalar. To proceed further we now 
assume that the fluid flow remains geodesic (i.e. A a — 0) despite the magnetic presence. In other words, 
we impose the force-free condition on the magnetic field (i.e. e abc _B b curLB c = = ^D a PB + divll a ). 

Appendix B. Calculating the limits 

Our method provides a small refinement of [11], allowing us to get slightly stronger limits (by about a 
factor of about two). In [11], limits on a ab were calculated in the following way: First, Eq. (A. 7) is solved 
for a ab , and then the separate limits in Eqs.(7)-(22) of [11] were inserted to give the following limit: 

J^l < |e a + ^ + Set + |ej. (R1) 

However, if, after solving Eq. (A. 7) for a ab , we use Eq. (41) of [8] to convert n ab — > r ab etc. (i.e., 
r a ~ iq a /Ap, r ab ~ 15n ab /8fi...), and then use Eqs. (l)-(4) of [11] (after expanding all derivatives, and 
using the relevant evolution equations), some terms cancel, and we get the tighter limit: 

M<4 + b4 + |4. (b.2) 

This gives simpler limits than using the method in [11]. As a further example, consider the limits on E ab 
in the case of no magnetic field (including cosmological constant): 

\E a b\ 



< 



( 18 27 1 4 

H 1 10el + 15e{* + 2e* + 3e** + — e\ + — e\* j + — HQ R e 2 , (B.3) 



e 

which is less than the corresponding limit in Eq. (28), [11]. 
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Our first problem is to find limits on H a b- In the absence of magnetic fields, limits on E ab may be found 
directly from Eq. (A. 9), using Eq. (A. 7). However, with magnetic fields present, this will give limits on 
the combination \E ab — ^H a b\, so we have to find separate equations for E a b and n a f,. We can do this by 
solving Eq. (A.13)and the time derivative of Eq. (A. 9), which gives 



IU= - — 



2Q t7ab 2 aab 



3 e 2 4 e 2 e 2 e 2 



+ + TTab + ^Q D (a%> - ^@ CUrl H ab, (B.4) 

with a similar equation for E ab . We have used Eq. (A. 4). In the linear regime, we can drop the 
angled brackets on time derivatives of PSTF tensors. Decoupling these quantities has introduced extra 
uncertainty into our equations, in the form of a ab . In Eq. (B.4), all the shear terms may be found in terms 
of the TA e simply by taking appropriate derivatives of Eq. (A. 7) and using Eqs. (11) and (12), followed 
by (5). The only term left to worry about is the curliJ ab term; however, we may use D c of Eq. (A. 15), 
given that 

\B c H ab \ < \B a B b a cd \ + \B a B b u c \. (B.5) 
Limits on the second gradient of the shear may be found from Eq. (A. 7): 

\B a B b a cd \ < ^H 3 {10561" + Ue» + 21$* + 27et ft } . (B.6) 
The rotation term may be found using u a = —curl B a p/2p = e abc B b B c p/8Hp and Eq. (A. 6): 

|D a D^ c | < E ff 3| 5c ttt + 5e ttt* + 6c tm}. (B.7) 
So we finally have 



\Jlab 

e 



25el + 5e* + -^e 



7~ 3 



t 8 6 . 18 + ' 
9,1 + ^ 2 + - b e 2 + y 4 



+ io e i + 7 P* + Pr + k + 



45 



15 



5e\ + e* 
9 



+ 



117 



tff 



27 



ttt* 



3e 



tt 



2 62 



\Egb 

e 



i 



i 



25 e t + 5e* + y e\ 



51 Jitt 
10 2 

9 t 



7 63 



18 

7 

si^ttt 

14 3 



4 + 



135 
— i 
14 



t* 



+ 



81 



(Bi 



+ iM + + 



5c 



135 



t* 



45 



t** 



27 



27 



243 
28"' 



t* 



81 



+ 28 63 



+ il! e ttt + ^Z e ttt* + 3 £ tt + 9 £ tt* + |i e tttt + |i e ttt 



8 -1 ' 2~ z ' 4~ z ' 20' z ' 28~ a 1 ,L> ' 
Comparing Eq. (B.9) with the limits found in the absence of a magnetic field, Eq. (B.3), reveals the extra 
complexity and uncertainty involved in having just one additional field. 



Appendix B.2. Inhomogeneity of the field 

To find limits on the inhomogeneity, we need to find limits on D a pe = — 3divII a . [Note that solving D a 
of Eqs. (A. 8) and (A. 16) for B a pB does not work since it does not separate gradients of pb and pm-] First 
we take D a of Eq. (B.4), and note that 

divcurl# a = icurl divi7 a = -^curl curl q a + i ( p M + ^p) curlw a , (B.10) 
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so 



|divcurlff a | < -|D a D bg | + - 



Pm + xA* |D Wb| 



(B.11) 



Hence we find that 



|Dq/0 B | 




< 



max |D a p B | 




} 



(B.12) 
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